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Abstract —This paper proposes a codebook for jointly 
quantizing channel direction information (CDI) of spatially 
correlated three-dimensional (3D) multi-input-multi-output 
(MIMO) channels. To reduce the dimension for quantizing 
the CDI of large antenna arrays, we introduce a special 
structure to the codewords by using Tucker decomposi¬ 
tion to exploit the unique features of 3D MIMO channels. 
Specifically, the codeword consists of four parts each with 
low dimension individually targeting at a different type of 
information: statistical CDIs in horizontal direction and in 
vertical direction, statistical power coupling, and instanta¬ 
neous CDI. The proposed codebook avoids the redundancy 
led by existing independent CDI quantization. Analytical 
results provide a sufficient condition on 3D MIMO channels 
to show that the proposed codebook can achieve the same 
quantization performance as the well-known rotated code¬ 
book applied to the global channel CDI, but with significant 
reduction in the required statistical channel information. 
Simulation results validate our analysis and demonstrate 
that the proposed joint CDI quantization provides substan¬ 
tial performance gain over independent CDI quantization. 


Index Terms —Three-dimensional (3D), multi-input-and- 
multi-output (MIMO), limited feedback, codebook design 


I. Introduction 

To meet the ever-growing data demand in future 
5th Generation (5G) cellular networks, one of the 
promising ways is to increase the number of antennas 
at the base station (BS) [|TJ. However, equipping 
large number of antennas at a BS is challenging due 
to the physical space limitation. This naturally calls 
for three-dimensional (3D) multi-input-and-multi- 
output (MIMO) systems [[2-51, where active antenna 
elements are placed in a two-dimensional (2D) array, 
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e.g., rectangular and circular arrays, or even non- 
planar arrays, at the BS. 

The non-linear arrays can provide both horizonal 
and vertical spatial resolution, and thus can support 
3D beamforming [|3j. 3D beamforming can be user- 
specific, which improves the signal to noise ratio 
(SNR) meanwhile generating less interference to ad¬ 
jacent users. Recently, a 3D MIMO prototype system 
operating at millimeter-wave bands was reported to 
support a multi-Gbps data rate service in macro cells, 
which provides high array gain to compensate the 
severe path loss [5J. 

In practice, the promised performance of 3D 
beamforming largely depends on how accurate the 
channel direction information (CDI) is obtained at 
the BS. In time division duplexing (TDD) systems, 
the CDI obtained by channel estimation in uplink 
can be used for beamforming in downlink if the 
antennas at the BS are perfectly calibrated 0. where 
the performance is limited by pilot contamination 
[6|. In frequency division duplexing (FDD) systems, 
limited feedback is widely used, where the CDI is 
firstly quantized at the user and then fed back to the 
BS |7J. Yet the feedback overhead is supposed to 
increase with the number of antennas [|8j, which is 
not acceptable for large antenna array systems. 

Spatial correlation is observed very typical in 
MIMO channels [9), due to small spacing between 
adjacent antennas and low angular spreads. In [T0|, 
B. Clerckx, el al. found that spatial correlation can 
be exploited to reduce the overhead for feeding back 
CDI significantly. The same conclusion was drawn 
in pTHl3j for spatially correlated massive MIMO 
channels, which indicates that FDD is also applicable 
for large antenna array systems without heavy feed¬ 
back overhead as supposed to be. 

Various codebooks have been proposed for spa¬ 
tially correlated channels. Theoretically, Llyod al¬ 
gorithm [14] can be applied to generate codebooks 
for 3D MIMO channels, but they are difficult to 
be used off-line in practice. There are some code- 
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books obtained from Llyod algorithms with reduced 


complexity, e.g., local-packing codebook [151, and 


gain and phase separated quantization [16]. A well- 
known codebook, rotated codebook [ 1/71, transforms 
the codewords optimized for uncorrelated channels 
(e.g., Grassmannian subspace packing (GSP) code¬ 
words) by channel correlation matrix. It was proved 
that the rotated codebook is asymptotically optimal 
in quantizing any spatially correlated channels as 


the codebook size becomes large [ 18[. An extension 
of rotated codebook to multiuser MIMO system is 


To design a desirable codebook in 3D MIMO 
channels under various angular spreads and reduce 
the amount of required spatial correlation informa¬ 
tion, we propose a joint CDI quantization strategy 
to quantize the CDI of two directions jointly. The 
designed codeword has a special structure led by 


the Tucker decomposition [25], which exploits the 


provided in 191, where not only each user’s own 
but also the other users’ correlation matrices are em¬ 
ployed for the codeword rotation. When designing 
codebooks for real-world cellular systems, practi¬ 
cal limitations such as constant modulus and finite 
alphabet need to be taken into account. Discrete 
Fourier transformation (DFT) codebook meets these 
limitations, which is suitable for highly correlated 


channels with uniform linear array at the BS [20]. For 
3D MIMO system with uniform rectangular array 
(URA), a Kronecker-product DFT codebook was 


proposed in [2T|]. However, as shown in [20], the 
performance of these pure DFT based codebooks de¬ 
grades severely when the angular spreads of channel 
increase. 

Generally speaking, there are two straightforward 
strategies to extend existing codebooks to 3D MIMO 
systems. The first strategy is global channel quan¬ 
tization, which expresses the 3D MIMO channel as 
a larger global vector, and reuse existing codebooks 
proposed for 2D MIMO channels. For example, the 


codebooks in [15 -17 22] can be applied to quantize 
the global 3D MIMO CDI directly. However, due to 
the high dimension of the resulting CDI vector, it 
is challenging for the BS to obtain accurate channel 
statistical information required by the global channel 
quantization to reuse existing codebooks, e.g., chan¬ 
nel correlation matrix used by the rotated codebook, 
which is known as the “curse of the dimension¬ 
ality” [23]. Moreover, large dimension incurs high 


unique features of 3D MIMO channels. 

The first feature of 3D MIMO channel is the inher¬ 
ent geometrical structure inherited from the regular¬ 
ity of antenna arrays. For example, the array response 
of 3D MIMO with URA for each ray in the channel 
can be decomposed into two subarray responses of 
lower channel dimensions in horizontal and vertical 
directions. Similar decompositions can be also found 
for those large arrays nested from several smaller 
identical subarrays. This feature can be exploited for 
coping with the high-dimension quantization prob¬ 
lem in large antenna arrays. 

The second feature is power coupling. When the 
array responses of 3D MIMO channels are decom¬ 
posed into two subarray responses respectively in 
horizontal and vertical directions, a common ray gain 
is shared by the two subarray responses and not de¬ 
composable. The power coupling is important in the 
codebook design, since it connects with the subarray 
responses from different lower channel dimensions. 
The feature can be exploited to avoid the redundancy 
in the CDI quantization. 

Although many research efforts have been made 
for the dimension reduction problem, they mainly 
employ the singular value decomposition (SVD) of 
correlation matrix to exploit the channel correlation, 
e.g., [26], As far as the authors known, the unique 


complexity in matrix operations [24]. The second 
strategy is independent CDI quantization, which 
quantizes the CDI of 3D MIMO in horizontal and 
vertical directions independently and reuses existing 
codebooks in each direction 0. This strategy is 
simple and avoids the problem caused by the high 
dimensional channel vector, which however comes at 
a cost of low quantization accuracy when the angular 
spread is large, as will be clarified later. 


features in 3D MIMO channels are not yet exploited. 
In this paper, we propose a codebook for jointly 
quantizing the CDIs in horizontal and vertical direc¬ 
tions, where a new codeword structure is introduced 
to reduce the dimension by using the features of 3D 
MIMO channels. The contributions are two-fold: 

1) A joint CDI quantization codebook is pro¬ 
posed to quantize spatially correlated 3D 
MIMO channels, which has a special code¬ 
word structure. The spatial correlation infor¬ 
mation required by the cobook for general 
antenna arrays is provided by using the Tucker 
decomposition. We show that the proposed 
codebook provides better quantization perfor¬ 
mance than the independent CDI quantization, 
and requires much less spatial correlation in- 
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Fig. 1. An example of 3D MIMO system with 8x8 URA at the BS. 


formation than the globally rotated codebook. 

2) The performance of joint CDI quantization 
codebook is analyzed. We show that the per¬ 
formance is the same as the globally rotated 
codebook for 3D MIMO systems with URA 
under the channels with weak identically in¬ 
dependent distributed (i.i.d.) rays. 

The rest of this paper is organized as follows. 
In Section II, we present the 3D MIMO channel 
model and channel representations. In Section III, 
we propose the joint CDI quantization codebook 
for spatially correlated channels. In Section IV, a 
solution to find the spatial correlation information 
for general antenna arrays is provided by using the 
Tucker decomposition. The performance of the joint 
CDI quantization codebook is analyzed in Section V. 
Simulation results are provided in Section VI and the 
paper is concluded in Section VII. 

Notations : (-) r , (•) /Y and (•)* are respectively the 
transpose, Hermitian and conjugate operation, 0 and 
0 are respectively the element-wise and Kronecker 
matrix product, | • [ and || • || are respectively the 
absolute value and norm, E{} means the expectation 
operation, diag(ir) is the diagonal matrix with diago¬ 
nal entries given in the vector x, and diag(X) is the 
diagonal matrix with the diagonal equal to that of the 
matrix X. 


II. System and Channel Models 

Consider a downlink 3D MIMO system where an 
array with N t antennas is mounted at the BS |2] 271, 
and all antennas are omni-directional. An example of 
3D MIMO system with URA is shown in Fig. |T] 

For easy exposition, we start by considering the 
3D multiple-input-and-single-output (MISO) system 


with single-antenna users, and then extend the de¬ 
sign to multi-antenna users. For simplicity, we will 
not distinguish 3D MISO channel from 3D MIMO 
channel hereafter. 

According to the generic spatial channel modeling 
methods |9j28]|29), the 3D MIMO channel consists 
of several scattering clusters distributed in the 3D 
space, where in each cluster there are multiple rays 
with small random angle offsets, which can be ex¬ 
pressed as a iVt-dimensional vector 


h = Yl , 1 9n, m a(Q 

L —'n=l L — *m=l 


(1) 


where g n>m e C is the random gain of the mth 
ray in the nth cluster with zero mean, the angle 
vector @ n , m = {9n t m, (f) n ,m} T is specified by the 
angle coordinates 9 nrn and o n , rn in the 3D space, 
and a(0„ im ) e C /VtX 1 is the corresponding array 
response. 

The expression of array response a( 0 n/m ) depends 
on the specific form of the array mounted at the BS. 
The array responses of some well-structured arrays 
can be decomposed into the Kronecker product of 
subarray responses, while the others may not. Usu¬ 
ally, the large antenna arrays nested from smaller 
identical subarrays are decomposable. For example, 
1, the array response of URA can 


as shown in [28 


be decomposed into two subarray responses respec¬ 
tively in horizontal and vertical directions as 




( 2 ) 


with 

a h (0n,m) = [l,e j2 ^ cosd "’ rn 
a v (</>n,m) = [ l,e j2vr ^ 


. . . (Vt-l)cos0 n , m jT 


? 1 

— IT O j^^\ cos 4>n,m . . . e j2n^(N v — l)cos<£„, m jT 


where d h and d v are respectively the antenna spacing 
in horizontal and vertical directions, N h and N v are 
the number of antennas at the URA in horizontal 
and vertical directions, A is the carrier wavelength, 
cos 9 n>rn and cos o nm are the direction cosines of 
the mth ray in the nth cluster respectively in the 
horizontal and vertical directions. 

One example of non-decomposable array response 
is the uniform concentric circular array (UCCA). 
Denote J and L respectively as the number of rings 
and the number of antennas equally placed on each 


ring in the array. As shown in [30], the array response 
of UCCA can be expressed as 


o(0n,m) = K</?l) , • • ■ Cl{Vl) } 


T IT 


(3) 
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with 

a(pf = [e j 2n ^f cos (<t>",m-‘Pi) cose n ,m 

g j 2-n- cos (</>„, m - If I ) cos e n , m j T 

where dj and ip t = 21"/L are respectively the radius 
of the jth ring and the /th radial direction in the 
array, j = 1,..., J, l = 1,..., L, 6 n , m and 0 n , m are 
respectively the directions of the mth ray in the //th 
cluster with respect to the positive x- and y-axis in 
the array. 

The CDI of 3D MIMO channel is h = h/\h\, 
which is of unit-norm f7j. In limited feedback MIMO 
systems, the CDI is quantized at the user by using a 
pre-determined codebook and then fed back to the 
BS. The quality of CDI available at the BS largely 
depends on the codebook. A desirable codebook 
should be judiciously designed by taking the channel 
features into account. To observe the unique features 
of 3D MIMO channels, besides the vector expression 
h e C NtXl , we also express the channel in matrix 
form. 

Taking the URA as an example, the 3D MIMO 
channel and the array response can be expressed as 

H G C NhXNv and A n>m = a h (e njm )a v ((j) nirn ) T (4) 

where h = vec (H), a( © n , m ) = vec(A ntm ), and 
vec(-) denotes the operation of vectorizing a matrix 
into a vector. Then, the CDI can be expressed as 
H = H/\\H\\. 

The matrix representation helps us to identify the 
different types of channel information in 3D MIMO 
channels with the URA, since the columns and rows 
of H respectively stand for the sub-responses at the 
horizontal and vertical directions. 

If the ray responses of an array are decomposable, 
i.e., o( O n ,m) = x Cg) y, we can easily find the matrix 
expression for each array response as xy 1 and the 
matrix expression for the channel, where x and y are 
not necessarily the sub-responses in the horizontal 
and vertical directions. However, the matrix expres¬ 
sions for arbitrary antenna arrays can not be obtained 
straightforwardly and will be discussed in Section IV. 

III. Codebook Design 

In this section, we propose a joint quantization 
codebook for spatially correlated 3D MIMO chan¬ 
nels, given that the channel matrix expression H is 
available. 


When tens or hundreds of antennas are placed 
in the array, say N t = 64 or 256, the dimension 
of vectors a(0 n m ) and h can be very large, and 
the dimension of corresponding channel correlation 
matrix is much larger since 

R = E{hh H } (5) 


which is of size N t x N t . In fact, the dimension 
N t x N t is large even for small value of Nf, and N v , 
say N t = N h N v = 8 2 for N h = N v = 8. Such a 
large dimension not only increases the computational 
complexity for MIMO signal processing but also 


results in the difficultly of CDI feedback [311. 

The channel correlation matrix R is important for 


many codebooks [17 19 26]. For example, when R 


is available, the rotated codebook can immediately 
be applied to quantize the 3D MIMO CDI vector, h. 
Specifically, the globally rotated codeword can be 
constructed from an instantaneous codeword g 0 of 

size N t x 1 as 03 


c g = R L/2 g 0 


( 6 ) 


where the codeword c g is of unit-norm. 

The globally rotated codebook has been proved 
to be asymptotical optimal to quantize arbitrary spa¬ 


tially correlated channels [ 181. Yet it is challenging 
to obtain the channel correlation matrix R at the BS 
since N t is large in 3D MIMO systems. In the sequel, 
we strive to reduce the dimension of quantization by 
exploiting unique features of 3D MIMO channels. 


A. Different Types of CDI of 3D MIMO Channels 

We begin with the following proposition to iden¬ 
tify different types of CDI in 3D MIMO channels. 

Proposition 1: Any 3D MIMO channel matrix H 
can be decomposed into H = UhH t U such that 
E{H t H?} = diag(A h ) and E{HjH* t } = diag(A„), 
where U h and U v are unitary matrices, X h and A,, are 
vectors with nonnegative entries. 

Proof: The proposition is easy to show by taking the 
SVD to channel matrix, which is given here simply 
for introducing notations. 

Denote the SVD of the left and right correlation 
matrices of 3D MIMO channel respectively as 

Rh = E{HH h } = U h dmg(\ h )U* (7) 
R v = E{H T H*} = £7„diag(A„)E7? (8) 
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where the entries of A/, and A,, are in a descending 
order. Considering that H t = U]] HU*„ and by using 
(|Tj) and ([8]), we can obtain the proposition. ■ 

This proposition suggests that for arbitrary 3D 
MIMO systems under arbitrary channels , we can 
always obtain two unitary matrices shown in ([7]) and 
([8]). Although this is nothing but taking the SVD to 
the 3D MIMO channel matrix, but the geometrical 
structure of the array has been explicitly exploited. 

Note that we can also obtain a unitary matrix U 
for the full channel correlation matrix by the SVD as 
R = Ud'\i\g(X)U n . To emphasize the difference, the 
unitary matrix U is referred to as statistical direction 
information, and the two unitary matrices U /,. U v are 
referred to as statistical sub-direction information. 

The statistical direction information can transform 
the 3D MIMO channel matrix H into at most NhN v 
uncorrelated channel gains expressed in the diagonal 
matrix diag(A). In contrast, the two statistical sub¬ 
direction information independently transform the 
channel matrix H into at most Nh and N v uncor¬ 
related channel gains respectively expressed in di¬ 
agonal matrices diag(Ah) and diag(A„). One desir¬ 
able property of statistical sub-direction information 
is that the dimensions of the correlation matrices 
of Nh x Nh and N v x N v are much smaller than 
NhN v x NhN v for statistical direction information. 
In the sequel, we refer to the two sub-directions as 
“horizontal” and “vertical” directions, although such 
notions only agree with their physical meanings for 
theURA. 

Denote j as the the average channel gain of 
the (i,j )th element in H t given by Proposition 1, 
e.g., Nfj = E{\H t)i f 2 }. By expressing u { and Vi 
respectively as the /'th column of matrix U /, and U v , 
then A t j can be also expressed by 

A. iJ =eH\u?Hv*\ 2 }. ( 9 ) 

which means A t J is the common average channel 
gain shared by the Ah statistical horizontal direction 
and jth statistical vertical direction. 

Different from X h in ([7]) and A, in ([8]) which are 
obtained separately, A is interacted jointly with U h 
and U v , and can be interpreted as the statistical power 
coupling information. 

B. The Proposed Codeword Structure 

Based on Proposition 1 and the observations in 
last subsection, we propose a joint quantization code¬ 
word structure for quantizing the CDI of 3D MIMO 


channel in matrix form, H = H /\\H\\, by separating 
the statistical sub-directions, the average and instan¬ 
taneous channel gains, which is 

Cj = U h {koG)U T v (10) 

where the codeword C j is normalized to have unit 
norm. By using vec(ABC T ) = (C <8> A)\ec(B), 
the codeword structure for quantizing the CDI of 3D 
MIMO channel in vector form h can be expressed as 

cj = ( U v <8> U h )dmg(X)g (11) 


where cj = vec(C'j), A = vec(A) and g = vec(G). 

The role of each part of the codeword is as follows, 

. unitary matrix U h is of size Ny l x r^, which tar¬ 
gets at the statistical sub-direction information 
Uh in horizontal direction (1 < r h < N h ), 

• unitary matrix U v is of size N v x r v , which tar¬ 
gets at the statistical sub-direction information 
U v in vertical direction (1 < r v < N v ), 

• nonnegative scaler matrix A is of size r/, x r v , 
which targets at the statistical power coupling 
information A between the horizontal and verti¬ 
cal directions, 

. instantaneous codeword G is of size ry, x r v , 
which quantizes the instantaneous power cou¬ 
pling information. In fact, the instantaneous 
power coupling information reflects the instan¬ 
taneous CDI of H t in Proposition 1. 


With the codeword structure in ( flO] ), the CDIs in 
horizontal and vertical directions are jointly quan¬ 
tized together with the power coupling. Therefore, 
we refer the codebook of codewords with this struc¬ 
ture as the joint CDI quantization codebook. 

Compared with the codeword structure of the 
globally rotated codebook shown in ([6]), the proposed 
codeword structure in <HD needs low-dimensional 
statistical channel information for rotation. 

In practice, the statistical information Uh, U v and 
A can be obtained at the BS either by uplink channel 
estimation or by feedback. 

In FDD systems where the downlink and the 
uplink are operated at separated frequency bands, 
estimating the downlink channel correlation matrix 
from uplink training symbols is still possible [ 321. 
However, such a method may not be used to estimate 
the channel correlation matrix R for 3D MIMO 
systems, which is of high dimension. The problem of 
estimating the correlation matrix for high-dimension 
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random vectors is recognized as “curse of dimen¬ 
sionality” and far from trivial p4| . This is because 
the dimension of the channel correlation matrix may 
be comparable with the number of collected training 
symbols such that the widely-used “sample covari¬ 
ance” estimation becomes invalid. 

The matrices Uh, U v and A can also be quantized 


and fed back [33J, 34], Since the spatial correlation 
information can be fed back in wide-band and in 
long-term, the feedback overhead can be almost ig¬ 
nored compared with the feedback for instantaneous 
CDI. This is especially true when the array feature 
is taken into account. For example, for the URA, 
according to Szego’s theory of Toeplitz matrices, 
the statistical sub-direction information U /, and U„ 
become a subset of DFT matrix when the size of the 
array increases GD- This will significantly simplify 
the feedback for the correlation information. 

Since our focus is to reduce the dimension of 
quantization by introducing new codeword structure, 
we assume that the statistical channel information 
U /,, U r and A are available at the BS. Then, the ma¬ 
trices used to construct the codewod with structure in 
m U h ,U v and A can be obtained as 

Uh. [til, • • • , Mr;,] 5 U v [tq , ■ • • , £V„] 

A ij = A ij, with 1 < i < r h , 1 < j < r v (12) 

The codeword quantizing the instantaneous CDI 
of 3D MIMO channel, i.e., G, can be obtained by 
using the codebooks designed for uncorrelated chan¬ 
nels with dimension of rh x r v . For example, GSP 
codebooks proposed in |7j can be used, which is 
optimized for uncorrelated Rayleigh channels. 

The dimensions r/, and r v in ( [TO] ) can be designed 
to further reduce the dimension of codeword by 
discarding trivial statistical sub-directions in strong 
correlated channel. We leave this topic for future 
study, since similar idea has already been developed 


in [26], and moreover, the dimension parameter can 
be optimized by considering various system param¬ 
eters, e.g., the error tolerance, storage space, and 
overall quantization performance. 


C. Extension to Multi-antenna Users 

When considering multiple receive antennas at the 
user, it is reasonable to assume that the statistical 
channel information in ( |T()| ) associated with differ¬ 
ent receive antennas are identical due to the small 
space separation. However, the instantaneous CDI 


associated with different receive antennas may differ. 
Based on this observation, we can easily extend the 
codebook design to the systems where each user is 
equipped with N r antennas. 

Denote the super-script in parentheses as the index 
of receiver antenna, the joint quantization codeword 
in CD for multi-antenna users can be given as 

(l,...,7Vr) _ r (1) (Nr), 

L J — l c J ) Ej \ 

= (U v ®U h )dmg(\)[g( 1 \---,g^] 

where g^\ ■ ■ ■ . g Nr) are respectively the instanta¬ 
neous codewords for different receive antennas. 


IV. Statistical Information for 
Arbitrary Antenna Array 

The proposed codeword structure is applicable for 
any antenna array if the channel can be expressed 
in matrix form as H. This is because Proposition 1 
is valid for arbitrary 3D MIMO channels. Although 
quite natural to the decomposable arrays such as the 
URA, the way to express the channels in matrix form 
for arbitrary antenna array is not obvious because of 
two issues. 

The first issue is the choice of dimensionality for 
the 3D MIMO channel matrix, i.e., Nh and N v with 
N h x N v = N t . Taking the UCCA array with N t = 32 
as an example, we may express the channel vector h 
in a matrix form H of size 8 x 4, 4 x 8, or the others, 
but we are not clear which expression achieves a 
better quantization performance because the array is 
not rectangular. 

The second issue is the “antenna grouping” in 
the channel matrix, i.e., which subset of antenna 
responses should be in the same row or column in 
H for a given dimensionality of N h and N v . The 
antenna grouping determines the values of Uh, U v 
and A, and eventually affects the performance of the 
proposed codebook. 

As a consequence, we need to solve these two 
issues when designing the joint codebook for 3D 
MIMO systems with arbitrary antenna array. The 
challenges for addressing the two issues are different. 
The choices satisfying N t = N h x N v are always 
limited, and thus exhaustive searching is efficient to 
find the best dimensionality. However, the possible 
choices of antenna grouping exponentially increases 
with N t , which belongs to a combinational prob¬ 
lem and is of prohibitive computational complexity. 
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Therefore, we focus on the second issue in the fol¬ 
lowing, i.e., to find proper U h , U v and A for a given 
N h and N v . 


A. Optimizing Uh, U v and A 

To find the desirable statistical information, be¬ 
fore solving the antenna grouping problem, we first 
reconsider the proposed codebook structure. Recall 
that the entries of H t defined in Proposition 1 reflect 
the instantaneous channel gains, and the GSP code¬ 
book can be used to quantize the instantaneous CDI, 
which is optimal for i.i.d. channels. Therefore, the 
proposed codebook will be optimal if the entries of 
H t are uncorrelated and the statistical information is 
perfect (i.e., JJ h — U h , U h — U v and A = A). 

We can reconstruct a full channel matrix R from 
the matrices with lower dimensions as 

R = E{vec(U h H t Ul)\ec{U h H t U^) H } 

= (U v © U h )E{vec{H t }vec{H t } H }(U v © U h ) H 

= (U V ®U ;i)diag(vec{A ©A}) (?/„©?/ h ) H (13) 


Problem ( [14] ) belongs to a classic approximation 
problem of Tucker decomposition [25], where X t is 
the core tensor, U h and U v are respectively the factor 
matrices with columns u, and v, as tensors. The main 
objective of Tucker decomposition is to decompose 
a higher dimensional matrix into low dimensional 
factor matrices, and the tensor core encompass all 
the possible interactions among the low dimensional 
tensors in the factor matrices. In other words, Tucker 
decomposition is to reduce the dimension in the large 
matrix by finding the structure properties. Moreover, 
Tucker decomposition is a generation of the matrix 
SVD with several desirable features, such as orthog¬ 
onality, decorrelation and computational tractability. 
However, in general the problem of finding the solu¬ 
tion of Tucker decomposition is NP-hard, and there 
are few efficient algorithms in use. 

Note that although the problem in ( fl4| ) contains 
three variables, it can be simplified by only optimiz¬ 
ing Uh and U v without losing the optimality. This 
is because given U h and U v , the objective function 
satisfies 


where ( [T3| ) is achieved when the entries of H, are 
uncorrelated. 

If R = R, the proposed codebook is identical to 
the globally rotated codebook, since the codeword 

in ( [TT| ) becomes c,j = R 2 g and the same as in ([6]). 
Unfortunately, in general cases, R f R, and the 
proposed codebook becomes inferior to the globally 
rotated codebook, which is asymptotically optimal. 
To provide good quantization performance, it is rea¬ 
sonable to find Uh, U v and A such that R is as close 
to R as possible. 

For simplicity, we define A t = vec{A © A}. The 
problem of finding the optimal U h , U v and A, for 
arbitrary antenna array can be modeled as 


min \\R-(U v ®U h )diag(X t )(U v ®U h ) H \\ F 

At 

= min||(EA,® U h ) H R(U v ®U h )-diag(X t )\\ F (15) 

At 

= min||diag ((L/j,© U h ) H R{U V ®U h ))-diag(X t )\\ F 

At 

+ 11off ((U v © U h ) H R(U V ®U h )) ||| 

= 11off((£/ V ®U h ) H R{U V ®U h )) HI (16) 


where ( fl5| ) is due to fact that the norm is unitarily 
invariant, off(X) is the operation to the matrix X 
with all zeros on the diagonal, ( fT6] ) is achieved by 
the optimal A t for a given U h and U v , which is 


diag(A t ) = diag [{U v © U h ) ]1 R{U v © U h )) (17) 


min \\R-(U v ®U h )diag(X t )(U v ®U h ) H \\ 2 F (14) 

U /, ,U .At 

s.t. U*U v =I Nv , 

U»U h = I Nh , 

Xt y 0 

where I n is the identity matrix of size n x n, and 
x >- 0 means each element in x is larger than 0. 
The problem in ( |T4l ) does not directly optimize the 
antenna grouping in the channel matrix, which is 
unnecessary since when given the optimal solution 
of Uh, U v and A,, we can obtain the matrices used to 
construct the joint quantization codewods. 


It can be verified that the operation for computing X t 
in ( fTT| ) is identical to that for computing A with Q. 

Since ( [76] ) does not depend on the parameter A,, 
we only need to find U /, and U„ from a new opti¬ 
mization problem by replacing the objective function 
in < fl4] ) with ( fl6] ). With the optimized U /, and U v , we 
can immediately obtaining optimal A f by using ( fTT] ). 

B. Closed-form Solution for U /,. and U v 

Since a closed-form solution is more desirable for 
practical use, we consider a modified problem to 
find U h and U v . Specifically, we impose an extra 
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constraint on X t into the new optimization problem, 
which is given by 


X+ A, 


(18) 


where A„ >~ 0 and A/, >- 0. Then, the term inside the 
objective function of problem ( p~4| ) becomes 

(U V ®U h )dmg(X t )(U V ®U h ) H =B®C (19) 

where B = U v diag(X v )U^ andC = Uhdmg{Xh)U^ 
are positive semi-definite matrices. 

With the constraint in ( [J~8] ) and considering ( [T9] ), 
the optimization problem to find U /, and U v becomes 

\R - B ®C\\% ( 20 ) 


nun 

B,C 

s.t. 


N v xN v 


, and C E S NhXNh 


uct decomposition [35]. There are many other bene¬ 
fits to consider the Kronecker product decomposition 
here. By using such an decomposition, many struc¬ 
ture properties of R, such as symmetry, definite, and 
permutations, can be inherited by the matrices B and 
C. It is worthy to note that these structure properties 
are usually led by the regularity of antenna array, 
e.g., symmetry, and nested subarrays. 

The solution to the problem in ( [20] ) is given in [ 35 ]. 
To obtain the solution, according to [35], we need to 


/ -Ri,i 

Ri,n v \ 

*= i 

(21) 

V Rn v , i • • • 

Rn v ,n v ) 

where the (i, j)th block denoted as R, , is of size 
N h x N h . Then a rearranged matrix is generated by 

R= [vec(i?i,i),vec(i? 2 ,i),- • 

■ ,\ec(R Nv , Nv )] T (22) 


which is of size N 2 x Nfi. 


Denote the largest singular value, the correspond¬ 
ing left and right eigenvectors of the SVD to R 
respectively as a 2 , u and v. Then, as shown in [35], 
the matrices B and C are obtained as 


vec(B) = au , and vec(C) = av 


where § nxn is the space of positive semi-definite 
matrices with the dimensionality of n x n. 


The new problem in ( ]20| ) offers a suboptimal so¬ 
lution of Ui, and U v with closed-form for the Tucker 
decomposition, and is known as the Kronecker prod¬ 


rearrange the matrix R. Specifically, the matrix R is 
divided into N v x N v blocks as 


(23) 


Moreover, Since R is symmetric and positive semi- 
definite, B and C are also symmetric and positive 


semi-definite [35]. Thus, the matrices B and C can 
be regarded as the correlation matrices. 

By letting R v = B. R h = C, using the SVD, we 
can obtain U v and Uh immediately. Then, by using 
( fT7| ), we can obtain the optimal A t or A. 

It is easy to validate that the statistical direction 
information U is approximated by the lower dimen¬ 
sion matrices as U v ® U /,, and the average channel 
gains in A are approximated by X t . By introduc¬ 
ing the structure of joint quantization codebook, we 
employ a new correlation matrix constructed with 
reduced dimension that approximates the full corre¬ 
lation matrix as close as possible. Therefore, with the 
obtained correlation matrices, we can not only solve 
the dimension problem for rotated codebook in [ 17], 
but also for the codebook in [ JT9| , both rely on the 
channel correlation matrices. 

V. Performance Analysis 

Since the globally rotated codebook is asymptot¬ 
ical optimal and can serve as a performance upper 
bound for other codebooks, we analyze the per¬ 
formance of the proposed joint CDI quantization 
codebook by comparing with the globally rotated 
codebook in this section. 

As addressed in Section IV.A, if the entries of H t 
are uncorrelated and the statistical information used 
in the codeword with structure in ( |T0| ) is perfect, the 
joint CDI quantization codebook will perform the 
same as the globally rotated codebook. However, it 
is unclear whether the entries of H t are uncorrelated 
or not. In what follows, we show that for 3D MIMO 
systems with URA such an uncorrelated property is 
valid under very general channel conditions. 

To facilitate analysis and gain useful insights, we 
consider weak i.i.d. rays in 3D MIMO channels, 
which assume, 

1) the angles 9 n>rn are i.i.d. for V n and m, 

2) the angles are i.i.d. for V n and m, 

3) the gains g njn are uncorrelated for V n and m, 
where the word “weak” comes from the fact that 
the third condition only requires the gains being 
uncorrelated but not being identically distributed. 

The above assumption does not weaken our per¬ 
formance analysis for realistic scenarios. First of 
all, the assumption of uncorrelated ray gains follows 
from the commonly-used uncorrelated scattering as¬ 
sumption for multi-path fading channels, which is 
validated by many channel models [9j[28j[29j. Sec¬ 
ond, i.i.d. angles are observed very typical and thus 
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are adopted in the 3GPP channel models [29]]. Fi¬ 
nally, no specific probability density distribution as¬ 
sumption is imposed on the gains and angles. For ex¬ 
ample, the angles can be uniformly distributed [28], 
Gaussian distributed [9|, or log-normal distributed 


[29], 

In the following, we show that the proposed joint 
CDI quantization codebook achieves the same per¬ 
formance as the globally rotated codebook when the 
3D MIMO channels have weak i.i.d. rays and the BS 
is equipped with URA. 


Lemma 1: For a 3D MIMO system with URA, the 
channel with weak i.i.d. rays can be expressed as 
H = HU T V , where U v is a unitary matrix and the 
columns of H are uncorrelated with each other. 

Proof: See Appendix [A] ■ 

Lemma 2: For a 3D MIMO system with URA, the 
channel with weak i.i.d. rays can be expressed as 
H = Uf L H, where U /, is a unitary matrix and the 
rows of H are uncorrelated with each other. 

Proof: The proof is similar to Lemma !.■ 


Lemma 3: For a 3D MIMO system with URA 
under the channel with weak i.i.d. rays, the entries 
of H t in the channel decomposition H = UhH t Uj ; 
given by Proposition 1 are uncorrelated. 

Proof: See Appendix [B] ■ 

By using these lemmas for the 3D MIMO system 
with URA under channels with weak i.i.d. rays, we 
can immediately find the relationship between the 
joint CDI quantization codebook and the globally 
rotated codebook. 


Theorem 1: For a 3D MIMO system with URA 
under the channel with weak i.i.d. rays, the proposed 
joint quantization codebook with perfect correlation 
matrices performs the same as the globally rotated 
codebook. 

Proof: See Appendix[C]H 

It follows that for 3D MIMO systems with URA 
under the channels with i.i.d. rays, the proposed 
joint CDI quantization codebook is asymptotically 
optimal as the codebook size increases, which is 
promised by the globally rotated codebook | [T%) . 

It is worthy to note that the Theorem 1 provides 
a sufficient condition on the 3D MIMO channels 
where the joint CDI quantization and globally rotated 
codebook are identical. The entries in H t may not 
be uncorrelated when other antenna arrays ( for 


example, UCCA) are used. This means for general 
conditions, the two codebooks are not identical. In 
these cases, the joint CDI quantization codebook, 
though still applicable, becomes inferior to the glob¬ 
ally rotated codebook. 


A. Further Comparison With Globally Rotated 
Codebook 

As shown in 0. the globally rotated codebook 
employs the channel correlation matrix, which is 
of size NhN v x NhN v . By contrast, the joint CDI 
quantization codebook employs three lower dimen¬ 
sion correlation matrices, each respectively of size 
N h x N h for R h given by ([7]), of size N v x N v for 
R, given by ([8]), and of size N h x N v for A given 


by ( fl2] ). Apparently, a high dimension matrix leads 
to a high complexity in generating the codewords to 
quantize the statistical information or transforming 
the instantaneous codeword. By using the joint CDI 
quantization, the computational complexity is signif¬ 
icantly reduced. 

In addition, different size of correlation matrix 
needs different amount of information to estimate or 
feedback. Taking an 8 x 8 3D MIMO channel as an 
example, the globally rotated codebook requires a 
correlation matrix with 64 2 = 4096 entries, while 
the joint CDI quantization requires the correlation 
matrix only with 64 x 3 = 192 entries. Even when 
the correlation matrices R. Rh, R, are Hermitian 
and Toeplitz, the globally rotated codebook needs 63 
complex elements in R, while the joint quantization 
only needs 7 complex elements in Rh, 7 in R v , and 
i'h x r v scaler elements in A. This implies that the 
proposed joint CDI quantization book requires much 
less channel correlation information to be estimated 
or fed back than the globally rotated codebook. 


B. Comparison With Independent CDI Quantization 

Another strategy, independent CDI quantization, 
can be applied to 3D MIMO systems, where the 
CDIs in horizontal and vertical directions are quan¬ 
tized independently |4]. By nature, joint quantization 
is not a simple aggregation of independent quantiza¬ 
tion in two directions. 

Denote the codewords independently designed for 
horizontal and vertical directions as C//, and Cj v , 
respectively. In Q, the CDI of 3D MIMO is recon¬ 
structed by two independently designed codewords 
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as Cj = CihcJ v . Taking independently rotated code¬ 
book as an example, c Jh and c Iv are respectively 
obtained by normalizing R x J 2 Ch and R}J 2 c v , where 
Ch and c, are instantaneous codewords for horizontal 
and vertical directions. The independent CDI quan¬ 
tization is easy to implement and is compatible to 
existing 2D MIMO systems, but at a cost of quanti¬ 
zation performance loss. 

First, independent CDI quantization will not be 
optimal to quantize the 3D MIMO channel if the rank 
of channel matrix H exceeds one. The rank of H, 
denoted as rank(iT), has a close relationship with the 
scattering environment. If there are multiple well- 
separated clusters or rays in the channel, it is easy 
to see that rank(iT) is far larger than one. However, 
the independent CDI quantization codebook always 
yields rank(CY) = 1, which fails to match the rank of 
generic 3D MIMO channel matrix. This implies that 
independent CDI quantization is only appropriate for 
quantizing the rank one 3D MIMO channels with 
a single ray. By contrast, the row rank and column 
rank of the joint CDI quantization codeword C j are 
respectively r h and r v . By selecting proper values 
of r h and r v , we can design joint CDI quantization 
codeword for any 3D MIMO channel matrix H. In 
practice, the values of r h and r v may be selected 
less than the rank of H for the purpose of dimension 
reduction. In these cases, the performance loss can 
be minimized by reserving the statistical directions 
in f fTTj] ) with nontrivial average channel gains. 

Second, independent CDI quantization results in 
redundant quantization on the channel information, 
since the codebooks are separately designed for hor¬ 
izontal and vertical directions from their own per¬ 
spective. To see this, we can find that the power 
profile of uncorrelated clusters seen at horizontal 
direction and at vertical direction in the transformed 
3D MIMO channel after SVD are respectively the 
eigenvalues of R h and R v given by (|T|) and ([8]). The 
two power profiles are generally not independent 
with each other. However, with independent CDI 
quantization the two power profiles are quantized 
separately, which leads to a redundancy in the CDI 
quantization. By contrast, when using joint CDI 
quantization, such a redundancy can be completely 
removed by quantizing the power coupling informa¬ 
tion. As a consequence, joint CDI quantization is 
more efficient to quantize the 3D MIMO channels 
than independent CDI quantization. 


VI. Simulation Results 

In this section, we compare the performance of 
joint CDI quantization with existing codebooks for 
3D MIMO systems with planar array by simulations. 

We consider multi-user MIMO transmission, 
where K single-antenna users are served by the BS at 
the same time-frequency resource with zero-forcing 


beamforming [36]. In the simulation, the users are 
randomly selected, which is equivalent to be se¬ 
lected by Round-Robin scheduling. The users are 
with homogeneous SNR but with different azimuth 
and elevation angles. In order to put the results of 
different antenna array size within one figure, the X- 
axis is set as the receive SNR. For transmit SNR, 
an extra log 2 N t dB should be considered due to the 
array gain. The gains of each ray are i.i.d. complex 
Gaussian distributed with zero mean. The sum rates 
are computed by the Shannon formula by averaging 
over 10 3 channel realizations. 

In the following, we consider two different angular 
spread modelings for the channel given by (jT]). One 
is simplified, which allows us to run simulations 
for large scale system efficiently. The other is more 


realistic, which is given in [29] and allows us to 
obtain more reliable results. It should be noted that 
our conclusions hold for both channel models. 

In the simulations, all the instantaneous CDI code¬ 
words are generated using random vector quantiza¬ 
tion (RVQ), which is easy to generate while with 
performance close to the GSP codebook [j37j. Specif¬ 
ically, two 5-bit RVQ codebooks are used for quan¬ 
tizing the instantaneous horizonal and vertical chan¬ 
nel directions in the independently rotated codebook. 
A 25-bit RVQ codebook are used for quantizing the 
instantaneous channel direction information respec¬ 
tively in the joint CDI quantization codebooks, and 
in the globally rotated codebook. 


A. Results in Simplified Scenarios 

In the simplified angular spread model, the az¬ 
imuth and elevation angles in (JT]) are, 

9n,m — $0 + $n + (24) 

fin,m 00 + fin T $fi n ,m (25) 

where 6 0 = U(— 60°, 60°) and fi 0 = U(— 45°, 45°) 

are respectively the mean of azimuth and elevation 
clusters, and U (a, b ) denotes uniform distribution 
with range from a to b, 9 n and fi n are the the azimuth 
and elevation deviation of the nth cluster from their 
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mean, S9 n . m and 3o lum are the random offsets. We 
set iV = 12 and M — 20, and model the 9 n and o n 
as i.i.d. Gaussian distributed with zero mean and a 
variance of cr 2 , where a is a parameter representing 
the angular spread. We model S9 n>m and do nm as 
Lapalacian distributed with a root mean square of 1°. 
As shown by [9j, such a simplified channel model 
well preserves the spatial features of the channels 
modelled in [129 ]. 

Both low angular spread (a = 5°) and large 
angular spread (cr = 20°) scenarios are evaluated. 

We evaluate the average sum rate of a 3D MIMO 
system with the configuration N h and N v using the 
joint CDI quantization codebook (JQC). For compar¬ 
ison, the performance of globally rotated codebook 
(GRC) is provided, which quantizes a N t = N h x N v 
3D MIMO CDI vector as a whole using the full 
correlation matrix R in ([6]). The performance of in¬ 
dependent quantization with rotated codebook (IQC) 
is given, which quantizes an Nh x 1 CDI vector 
for each direction individually using the statistical 
information R h in ([7]) and R v in ([8]). 

1) Performance for 3D MIMO Systems with URA: 
We first evaluate the performance of joint CDI quan¬ 
tization codebook for 3D MIMO system using the 
URA, where N h — N v . 

The average sum rates of the 3D MIMO systems 
using different codebooks are shown in Fig. [2] where 
all the channel statistical information are perfect. It 
is shown that under different scenarios, the curves 
for the joint CDI quantization codebook overlap 
with those of the globally rotated codebook, which 
validates Theorem 1. Under the same angular spread, 
when the number of users grows, the performance 
gap between the joint and independent CDI quan¬ 
tization increases. Compared with the low angular 
spread, the performance gap is larger in the channel 
with large angular spread. 

The average sum rates of 3D MIMO systems using 
the URAs with different sizes are shown in Fig. [3} To 
show the impact of imperfect statistical information, 
we also evaluate the performance of the joint CDI 
quantization codebook with quantized statistical in¬ 
formation (labeled as “JQ w/DFT”). The statistical 
sub-direction information in horizontal and vertical 
directions are respectively quantized by a 8-bit DFT 
codebook j20j, where each codeword consists of r h 
and r v adjacent columns in the DFT matrix. The 
values of r h and r v are selected to quantize the first r h 



Fig. 2. Average sum rates of 3D MIMO systems using URA for 
different angular spreads and numbers of users, where Nh = N v = 8 
and B = 8. 



Fig. 3. Average sum rates of 3D MIMO systems using URA with 
different sizes, where K = 4, a = 20°, and B — 8. 


and r 
and 


y ^ r h ^ 

dominant directions such that h,t >0.9 

S^ rv \ ■ 

> 0.9, where A h,i and A„,i are the /1h el- 


sr^N v x 

Z^i= l 

ement in X h and A r respectively. The power coupling 
matrix are quantized by a 8-bit RVQ codebook. 

As shown in the figure, when the array size 
increases, the performance gaps of the joint CDI 
quantization codebooks with perfect and quantized 
statistical information keep almost constant under 
different SNRs and numbers of antennas. This is 
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Fig. 4. Average sum rates of 3D MIMO systems using UCCAs with 
different angle spread and number of bits, where K = 4 and a = 5°. 


because the DFT codewords is a natural choice for 
representing the statistical directions for 3D MIMO 
with URA, as shown in CCD- Moreover, when the 
array size increases, the performance gap between 
the joint and independent CDI quantization code¬ 
book increases. This is because the rank of channel 
matrix increases with the size of antenna array even 
for a fixed angular spread. It is seen in the figure, 
as the size of antenna array increases, the sum rate 
gap between perfect CDI and the globally rotated 
codebook are reduced. This is because the inter¬ 
user interferences are reduced in the large antenna 
array when the full spatial correlation information 
is available. It indicates the importance of exploiting 
the channel correlation information in the codebook 
design in large array systems. 

2) Performance for 3D MIMO Systems with 
UCCA: We then evaluate the performance of joint 
CDI quantization codebook for 3D MIMO systems 
using the UCCA at the BS. The configuration of 
UCCA is given by J = 8, L = 8, d\ = 0.5A, 
and dj = dj -1 + 0.5A for 1 < j < J. The 
average sum rates of 3D MIMO systems are shown 
in Fig. |4} When obtaining the statistical information 
by using the proposed explicit solution, we set the 
dimensionality as N h = N v = 8, and similar results 
can be observed by other settings. 

As shown in the figure, under different angular 
spreads and numbers of bits, the performance of 
joint CDI quantization codebook is inferior to that 


of rotated codebook but with acceptable performance 
loss. Moreover, the joint CDI quantization outper¬ 
forms the independent quantization significantly. 

B. Results in More Realistic Scenarios 

In this subsection, we consider more realistic sce¬ 
narios in the simulations. The azimuth and elevation 
angles in the channel are modeled by log-normal 


distributions as in [29], where two scenarios are 


evaluated: 3D urban micro (UMi) with non line-of- 
sight (NLOS), which is referred as 3D UMi, and 3D 
urban macro (UMa) with NLOS from an outdoor 
BS to an indoor user, which is referred as 3D UMa. 
The antenna spacing is 0.5A for both horizontal 
and vertical directions. Unless otherwise specified, 
the main 3D MIMO channel parameters are listed 
in Table I, where d is the distance from the BS 
to the user, h is the height of each user. We set 
d = 100 m for 3D UMi, and d = 250 m, h = 1.5 
m for 3D UMa. This channel modeling considers 
more practical issues, such as the distance-dependent 
channel statistics, and the coupling effect between 
the statistics like delay spread and angular spread, 
which is more complicated than the simplified one in 
the previous subsection. 

We evaluate the average sum rate of an 8 x 8 
3D MIMO system using the joint CDI quantization 
codebook with different codeword dimension: JQC 
with full dimension (r h = N h and r v = N v ); JQC 
with low dimension (r h = 2 and r v = 2). 

The average sum rates of the 3D MIMO systems 
with the URA are shown in Fig. [5] where 3D UMi 
and 3D UMa channels are respectively used. It is 
shown that the joint CDI quantization codebook with 
full dimension achieves the same performance as the 
globally rotated codebook. Even when using the joint 
CDI quantization codebook with reduced quantiza¬ 
tion dimension, there is still substantial performance 
gain over the independent CDI quantization code¬ 
book. In both scenarios, when the number of users 
K grows, the performance gap between the joint 
and independent CDI quantization increases. The 
results are similar to the ones obtained in previous 
subsection. 

Compared with 3D UMa, the performance gap 
between the joint and independent CDI quantization 
is larger in 3D UMi under the same number of users. 
This is because the scattering environment in 3D 
UMi is richer than in 3D UMa due to large variance 
given by the parameters of log-normal distribution. 










TABLE I 

Main Parameters of 3D MIMO Channels 
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Scenario 

3D UMi 

3D UMa 

Mean of azimuth clusters 

U (—60°, 60°) 

U{- 60°, 60°) 

Mean of elevation clusters 

U{- 45°, 45°) 

U(- 45°, 45°) 

Mean of log delay spread (DS) (log 10 ([s])) 

-6.89 

-6.62 

Variance of log DS (log 10 ([s])) 

0.54 

0.32 

Mean of log azimuth spread (AS) (log 10 ([°])) 

1.41 

1.25 

Variance of log AS (log 10 ([°])) 

0.17 

0.42 

Mean of log elevation spread (ES) (log 10 ([°])) 

max[—0.5, -2.1(d/1000) + 0.9] 

max[—0.5, -2.1(d/1000) - 0.01(/i - 1.5) + 0.9] 

Variance of log ES (log 10 ([°])) 

0.6 

0.49 

Number of clusters 

19 

12 

Number of rays per cluster 

20 

20 




(a) 3D UMi (b) 3D UMa 

Fig. 5. Average sum rates of 8 x 8 3D MIMO systems with B = 4. 


VII. Conclusions 


independent CDI quantization. 


We proposed a joint CDI quantization codebook 
for 3D MIMO systems in this paper. By exploit¬ 
ing the unique features of 3D MIMO channels, the 
proposed codebook is composed of codewords with 
a special structure, reflecting different types of in¬ 
formation in the 3D MIMO channel direction. The 
proposed codebook is easily adapted to different 
spatially correlated channels and is applicable to 
different forms of arrays. Analytical analysis indi¬ 
cated that the proposed codebook achieves the same 
performance as the globally rotated codebook for 
uniform rectangular array under very general chan¬ 
nel conditions but needs much less spatial correlation 
information. Simulation results validated the anal¬ 
ysis and shown substantial performance gain over 


Appendix A 
Proof of Lemma 1 


For notational simplicity, we express the Ah col¬ 
umn of 3D MIMO channel matrix H in ([Tj) and the 
jth column of its transpose H T as, 


kN 


yM 


Ci = 


Z— ( m= l 

M 


* 71=1 ■ 

.N 


= 


‘ ■ J n=l < J rn— 1 


Qn,mQ , v,i (0n,m )a h (0 n,m) 

9n,m a h,j{9 n,m )Cb v (0 n,m ) (26) 


where a Vji ((j) n , m ) and a h .j(O n:m ) are respectively the 
Ah entry of vector a v (6 nrn ) and the Ah entry of 
a h (9 n , m ),l<i<N h ,l<j<N v . 

We can set H = U h H t , where U h and H, are 
defined in Proposition 1. 
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Denote the zth column of H as c*, the A;th row of H 
as , and the (th entry of r k as r k Then, we have 


and 


E {'rk/tj} = 0, Vfc ^ £ and i ^ j, (34) 


r k = U"r k 


(27) 


The cross-correlation matrix between the zth and 
jth column of H for i / j can be obtained as 


because the off-diagonal elements of matrix A]f are 
zero. 

From ((281) and (J34b, we obtain 


(28) 


E{c,cf] = E {{r, 

For a 3D MIMO system with URA, for any 1 < 


k,P < N v and k ^ £, by using ( J26[ ) we have 

N M 

E { r krf} = E (EE*. 


n= 1 m=1 
N M 


EE 


E{cicf} = 0, Vz 7 ^ j (35) 

which means that the columns of H are uncorrelated 
with each other. The lemma is proved. 

Appendix B 
Proof of Lemma 3 

Denote the nth column of H t as c tn . By letting 
H t = U k H, i.e., c tn = U k c n , forn^m we find 
that 


n= 1 m= 1 
N M 


E{ c tnCu n } = U k E{c n cf n }U h = 0, Vn ^ m (36) 


tH i 


*H 


: E E^( I 5 '"’” 1 1 2 } E i a h,k {0n,m)a* h/ (0 n>m ) 

n=lm=l 

O'! ;(0n,m)®v (^n,m)} (29) 

iV M 


where the last equality is from ( [35] ) in Lemma 1. 

Similarly, denote the nth row of H t as r[ n . By 
letting H t = HU* V , i.e., r tn = U^r n , for n ^ m 
we find that 


n=lm=l 


E { r tnrf m } = U“E{r n r*}U v = 0, Vn / m (37) 

where the last equality is obtained from Lemma 2. 

From ( |36[ ) and ([37]), the entries in H t are uncorre¬ 
lated. The lemma is proved. 


£{a„(0n,m)af {<j>n,m)} (30) 

which are derived under the assumption of weak i.i.d. 


rays: ( |29| ) is because the gains <j n . m are uncorrelated 
with each other and also independent with the angles 


Appendix C 
Proof of Theorem 

From Lemma 3, the entries of H t are uncorrelated. 
When the spatial correlation matrix of the 3D MIMO 
E {a v ((p n ,m)ay (On.rn )} — Ry = U V A^U ^, for Vn and channels is perfect, from ( [T3| ) we have 

R = R = E{vec{H)vec{H) H } 


d n m and 0n iTn ; and ( p0[ ) is because the angle 9 n , m is 
independent with 4) n ,m- 

For i.i.d. variables 4> n ,m, we have 


(31) 


where the last equality is the SVD decomposition of 
positive semi-definite Hermitian correlation matrix 
R° u , and is a positive diagonal matrix. 

Then, from ( [3T)[ ) and ( |37] ), we have 


= (U v © U / t )diag(vec{A © A})(U V © U, 


\H 


E{r k r f } = U v AyU 


UttH 

V 


(32) 


where 


N M 


Af = A° v 2 }E{a h A0n, m )a* h/ (6 n , m )} 


n=lm=l 


is diagonal. 

From ( [27] ) and ([32]), we have 


E{r k rf} = U*E{r k rf }U V = Af (33) 


Then, the rotation matrix used in the globally 
rotated codeword in ([ 6 ]) is 

R i/2 = (U v © C/ ft )diag(vec{A}) (38) 

which is the same as the rotation matrix used in 
the joint CDI quantization ( ]TT[ ) when the statistical 
information U h = U h ,U v = U v , and A = vec(A) 
are perfect. 

Therefore, when all the instantaneous codewords 
for the joint CDI quantization codebook and those 
for the globally rotated codebook are identical, the 
two codebooks achieve the same performance. The 
theorem is proved. 
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